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Abstract By using rather conservative estimates based| 
on the simplest polar cap model, we search the ATNF 
Pulsar Catalogue for strongly magnetized stars that 
could accelerate rclativistic protons up to the curva- 
ture pion production threshold. The best candidate 
turns out to be the 16 ms pulsar J0537-6910, but the 
corresponding characteristic parameter x = a/m p is 
yet too small to give origin to observable signals. We 
show that, for pulsars with period P w 1 ms, a surface 
polar magnetic field B 10 12 G is required in order 
to induce detectable curvature pion radiation from ac- 
celerated protons in the magnetosphere. Some other 
emission processes are also considered. 
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1 Introduction 



stars, magnetic fields, 



Processes related to the emission of pions by accelerated 
protons like, for instance, p + — > jp + 7r°, have been con- 



ter characterizing these phenomena is the dimension- 
less quantity x = o,/m p , where a and m p stand, 
respectively, for the proton proper acceleration and 
rest mass. (Unless otherwise stated, natural uni- 
ties are adopted through this work.) For x <C 1, 
one can employ a semiclassical approach where the 
proton is described by a classical current while the 
pion is considered fully quantized field. This 

corresponds to the so called no-recoil approximation, 
for which many analytical formulas are available. 
(For so me recent works on the subject, see, f or in - 
stan ce. dHerpav fe Patkosl l2008t iHerpav et al. 1 l2008l) 
and (iFregolente. Matsas fc Vanzellall2006 ; Fregolente fc Saafr 
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for the electromagnetic and gravitational cases, 
respectively.) In particular, we expect a strong sup- 
pression of the pion emission for x ~ * since such a 
process is well known to be inertially forbidden . For 
Y ^> 1, a full quantized treatment is mandato ry (jRitus 
19851 : Berezinskv. Dolgov fc Kachelriess 1995 ). How- 



ever, numerical evidences suggest that some appro- 
_p_riatc limits of certain results obtained with the no- 



ma 

sidered in the literature since the six ties dGinzburg fc Zharkc-vj|nu ap proximation can still be c onsidered as good 



JJ6i;|Zharkoj|l965t[Grnzburg fc Sirovatski||l965fc |Ritus| estimates |Tokuhisa fc Kaiinol Il999l) . For the case of 



1985). The present available observational data on 



strongly magnetized astrophysical objects could pro- 
vide an exciting scenario to test theoretical predic- 
tions involving such kind of inertially forbidden pro- 
cesses. The relevant observer-independent parame- 
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the usual photon synchrotron radiation, for instance , 
quantum effects were carefully considered (jErberl 1 19661 ) 
and the semiclassical approach was found to be accurate 
within a few percent in the limit ^> 1, Notice that, 
in the cases where the proton acceleration is caused by 
strong magnetic fields, the parameter x can be written 
in a similar way to the usual synchrotron radiation: 
X = jB/B CI , where 7 is the Lorentz factor, jB is the 
magnetic field in the instantaneous reference frame of 
the proton, and B CY = m^/e « 1.5 x 10 20 G is a crit- 
ical magnetic field strength, denoting the limit of the 
validity of the no-recoil approximation. 

Thus, we have basically two regimes where the de- 
tails of the pion emission by accelerated protons arc 
known: the no- recoil regime a -C m^, where m„ w 135 
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MeV stands for the ir° mass, and the ultra-accelerated 
one, a ^> m p sa 938 MeV. The intermediate accelera- 
tion range corresponding to m T < a < m p , which could 
encompass interesting astrophysical applications, is still 
open for debate. It is clear, however, that for such accel- 
eration values, the pion emission channel will certainly 
be excited. In fact, the ir° radiation process starts to 
become competitive with the usual photon synchrotron 
radiation only for y w m^lm v { Tokuhisa fc Kaiinol 



2 Particle acceleration in the polar cap model 

We use here the simplest polar cap model for pul- 
sars. (For a compreh ensive review on the subject, see 



Harding fc Lai We recall that, in the magnc- 



tosphere of strongly magnetized pulsars, one expects an 
electric field component parallel to the magnetic field, 
which could accelerate charged particles up to relativis- 
tic velocities along the magnetic field lines near the po- 
lar region of the star. Pulsars with polar caps for which 
JIB < 0, where f2 and B stand, respectively, for the ro- 
tation and the magnetic field vectors, could effectively 
accelerate protons away from the polar region surface. 
Jly 80 callcd inner acceleration region (|Harding fcLaj 

by th e weak channel in muons and neutrinos (jNakamura et Jfl jP' the accelerated particles are expected to radiate 



1999t iHerpav et al. Il2008t ). 



The dominant decay mode for ir° is th e electromag- 
netic channel, giving origin to two photons ljNakamura et al 
2010l) . From the observational point of view, charged 



20ld) . P rocesses like 



r + have also been awav their transverse energy so efficiently that they re- 

.v fc PatkoskOQSUFregoJnte. Matsa/fc V»laf ectively constrained to move alon S the curved 

magnetic field lines with constant energy, correspond- 



20061) . and they have essentially the same two regimes 
according to the value of the parameter y. We no- 
tice, however, that the semiclassical analysis of such 
processes in magnetic fields is rather tricky since the 
neutron will not be accelerated and, consequently, 
will follow a different trajectory from the proton, im- 
plying that the initial (p + ) and final (n) currents 
are qualitatively different, in sharp contrast to the 
gravitational casertFregplente. Matsas fc Vanzella 2006 ; 
Fregolente fc Saal 12008?!). On the other hand, cases in- 



with m, 
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volving vector mesons like p + 

MeV, as those ones considered in Tokuhisa fc Kaiinol 
(1999). are analogous to the neutral n° case. 

In the present work, we searc h the ATNF Pulsar 
Catalogue ljManchester et al.ll2004j ) for strongly magne- 
tized pulsars that could give rise to appreciable values 
for the parameter y. We employ rather conservative 
estimates based on the simplest polar cap model for 
pulsars. The highest value of y we found corresponds 
to the 16 ms pulsar J0537-6910. Unfortunately, the ob- 
tained value y ~ 10~ 6 is too small in order to give ori- 
gin to observable signals. Nevertheless, our estimates 
reveal that a pulsar with period P = 1 ms and sur- 
face polar magnetic field B = 10 12 G would suffice to 
induce detectable pion radiation from accelerated pro- 
tons. Our results can be applied also to the weak inter- 
action channel emission, as in the process p + A ne + v. 
The threshold in this case correspond to the electron 
mass m e m 0.5 MeV. The y values for the known pul- 
sars, however, are not enough to excite this emission 
channels cither. For a P = 1 ms pulsar, a magnetic 
field B ?s 4 x 10 10 G would be required to induce the 
weak interaction emissions. 



ing to some Landau level. 

The curvature radius R c of th e polar magnetic field 
lines can be approximated by ([Harding fc Muslimov 
19981 ) 



R, 



8G 2 
9?r 



-rP, 



(1) 



where r is the radius of the star, P is the rotation pe- 
riod, and G is a curved space correction term, which 
tends to increase R c by 25% - 30% when compared to 
the flat space formula. We assume the flat space value 
(G = 1) since we are going to estimate the maximum 
possible centripetal acceleration attained by protons. 
The curvature radius ([1} does indeed correspond to the 
last open magnetic field line of the pulsar. For a rela- 
tivistic particle following a circular trajectory with ra- 
dius R c , the proper centripetal acceleration is given by 
a w 7 2 / R c - Assuming that the kinetic energy e p of the 
accelerated proton moving along the last open magnetic 
field line was obtained from the parallel electric field, 
one has 



7m p = rjeAcf), 



(2) 



where 77 ac counts to the efficiency of the acceleration 
mechanism) Zhang fc Harding 2000t Zhang et al. 2003 ) 
and 



A4> = 2^ 



(3) 



is the Goldreich - Julian maximum potent ial drop near 
the star surface (|Goldreich fc Julian! 1 1 9691 ) . with B be- 
ing the pulsar surface polar magnetic field. The 
maximum efficiency is expect to be as high as 0.85 
( Zhang fc Hardinj2000l) . but here we suppose total ef- 
ficiency, r) = 1. Finally, the maximum centripetal ac- 
celeration of such protons leaving the polar cap along 
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the last open magnetic field line can be approximated 

by 



(4) 



For our purposes here, it is convcnicntc to introduce 
the usual dimensionlcss quantities P_ 3 and B12 defined, 
respectively, by P = P_ 3 x lCT 3 s and B = B 12 x 10 12 G. 
Assuming a typical radius r = 10 4 m for the pulsar, one 
has finally 



= — « 0.35B? 2 PI 



(5) 



The total energy acquired by the proton is another im- 
portant quantity in our analysis. From ([2]), we have 



„ 2 B r3 
7 = 27T m„ 77 , 

leading that 

7 = 7 x 10 9 Bi 2 Pr 3 2 

for stars of the typical size. 

3 Results 



(6) 



(7) 



We evaluate x as given by (Et for all the 170 pul- 
sars in ATNF catalogue dManchester et alj 120041) with 
known P and B values. Table [1] shows the ten pulsars 
with the highest values of x- As one can see, these val- 
ues are too small to give rise to detectable signals, since 
the pion emission channel is effectively excited only for 
accelerations such that x ~ rn^/rrip ps 0.14. Never- 
theless, the proton can indeed reach considerable high 
energies, as the last column of Table [T] shows. For the 
J0537-6910 pulsar, for instance, the proton can reach 
an energy e p pa 2.3 x 10 5 GeV. We notice that for the 
pulsar with the highest surface magnetic field (J1808- 
2024, B w 2 x 10 15 G), one has only x ~ 5.2 x 10~ 12 , 
since its period is considerable large (P sa 7.6 s). 

Figure [T] show the (P_3,Pi2) pairs for the consid- 
ered pulsars. All the known pulsars are contained in 
the regions bounded by the curves ((5|) and |Zj) with, 
respectively, x — 10 -6 and 7 = 2 x 10 7 (upper) and 
X = 10~ 14 and 7 = 2 x 10 1 (bottom). The small val- 
ues obtained for x ar e, of course, a consequence of the 
large proton mass m p . Had we considered electrons in- 
stead of protons, we would have P cr pb 4.3 x 10 13 G, 
and the parameter x would be 10 orders of magnitude 
as larger. The regime x ~ 1 can be easily reached 
for electrons, explaining the rich physics of electrons 
in typical pulsar magnetospheres. Processes involving 



muons in circular orbits, like fj," — > e~P e v fi for instance, 
have also been investigated recen tly in the gravitational 
context (jFregolente k, Saal 120081 ) . The muon mass is 
rn p = 107 MeV, implying a muon acceleration param- 
eter x = a/m^ pb 7.4 x 10~ 4 for the pulsar J0537-6910. 
This is not enough to lead to any app reciable deviation 
of the muon half life as considered in iFregolente fc Saa 
( 20081) . 

Since we have x ^ lj we can a PPby safely the for- 
mulas obtained in the no-recoil approximation for the 
process p + A p + Tr°. The total energy emitted by the 
pion synchrotron radiation, or in other words, the pro- 
ton ener gy loss rate due to curvat ure pion radiation is 
given bv jTokuhisa fc Kaiino|[T999l ) 



9 2 ( VSnin 

Ac = —^m^mpxexp 

V3 \ X rn p 



(8) 



in the limit x 1, where g 2 pb 14 is the strong coupling 
constant. The so called cooling time t n = e p /lZ is an 
important parameter corresponding to the time scale 
for which the process becomes relevant. For our case, 
we have 



t 7T = 6x 10" 25 7x _1 exp (0.25x _1 ) s, 



(9) 



for stars of the typical size, with x an d 7 given, re- 
spectively, by ([5]) and ([?])• The exponential term in 
© gives rise to unreasonable high values of t v for the 
pulsars in the ATNP catalogue when co mpared, for in- 



stanc e, with the acceleration time t a (jHerpav et al 



2008| ). From Figure Q] we can also infer the location 
of the curves with constant t„ in the (P-3, P12) plane. 
Notice that the curves of constant x an d 7 are close in 
the log-log plot, since their inclination are, respectively 
9/4 and 2. Regions of approximately constant x and 7 
would give rise, according to ((9|) , to approximately con- 
stant t^. So, the curve of constant t n passing through 
the point (P_3,Pi2) will be quite close to the lines of 
constant x an d 7 intercepting at the same point. For 
the the values of x an d 7 corresponding to the upper 
bound of the known pulsars, we would have unreason- 
able highs values for t^. This is, of course, a conse- 
quence of the strong suppression of the pion radiation 
emitted by quasi-inertial protons. It is clear that no 
currently known pulsar could origin detectable signals 
of the curvature pion radiation emitted from acceler- 
ated protons in the magnetosphere. 
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Fig. 1 The values (P-s,B 12 ) for the 1700 objects in the ATNF Pulsar Catalogue l|Manchester et al.ll20oll ) with known 
period P = P_3 x s and surface magnetic field B = B12 x 10 12 G. The uppermost red (solid) and green (dashed) lines 
correspond, respectively, to the curves of constant \ w m^/mp (01 and 7 ~ 4.5 x 10 9 ([7]). The other red (solid) and green 
(dashed) lines correspond, respectively, to the lines with x ~ 10 -7 , 10 -9 , . . . , 10 -14 and with 7 ~ 10 7 , 10 6 , . . . , 10 2 . All the 
pulsars in the ATNF catalogue are in a region of the plane (P-3, -B12) from where no detectable signals are expected from 
curvature pion radiation emitted by accelerated protons in the magnetosphere. 



Pulsar 


P-3 


B\2 


X-7 


77 


J0537-6910 


16.12 


0.93 


11.0 


2.49 


B0531+21 


33.08 


3.78 


7.26 


2.42 


B0540-69 


50.50 


4.98 


1.88 


1.37 


J1813-1749 


44.70 


2.65 


0.92 


0.93 


J1400-6325 


31.18 


1.11 


0.82 


0.80 


J1747-2809 


52.15 


2.88 


0.54 


0.74 


J1833-1034 


61.87 


3.58 


0.39 


0.65 


J0205+6449 


65.59 


3.61 


0.30 


0.59 


J2229+6114 


51.62 


2.03 


0.28 


0.53 


J1617-5055 


69.36 


3.10 


0.17 


0.45 



Table 1 Top ten \ parameters, calculated accord- 
ing to B, for the 1700 object s of the ATNF Pulsar 
Catalogue (|Manchester et alj [2004) with known P and B 
values. The last column is the total proton energy given 
by ©. The pulsar period and magnetic field are, respec- 
tively, P = P-3 x 10" 3 s and B = B12 X 10 12 G. Analogously, 
we have x — X-7 x 10 -7 and 7 = 77 x 10 7 . 
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